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A method of and apparatus for processing seismic data 

The present invention relates to a method of processing multi-component seismic data. 
• It particularly relates to a method of processing seismic data to determine a calibration 
filter that calibrates one component of the seismic data relative to another component of 
the seismic data. The invention further relates to an apparatus for processing seismic 
data. 

Figure 1 is a schematic view of a seismic surveying arrangement. In this figure the 
smrveying airangement is a marine surveying arrangement in which seismic energy is 
emitted by a seismic source 1 that is suspended within a water colmnn 2 &om a towing 
vessel 3. When the seismic source 1 is actuated seisntiic energy is emitted downwards 
and is detected by an array of seiamic receivers 4 disposed on the seafloor 5, (As used 
herein the temi "seabed" denotes the earth's interior, and the term "seafloof denotes 
the surface of the seabed.) 

Many seismic surveys now use multi-component receivers that record two or more 
components of the seismic energy incident on the receiver. For example a 3-component 
(3-C) seismic receiver contains three orthogonal geophones and so can record the x-, y- 
and z-components of the particle motion at the receiver (the particle motion may be the 
particle displacement, particle velocity or particle acceleration or even, in principle, a 
higgler derivative of the particle displacement). In a marine seismic survey a 4- 
component (4-C) seismic receiver can alternatively be used. A 4-C receiver contains a 
pressure sensor such as a hydrophone in addition to three orthogonal geophones and so 
can record the pressure of the water colimm (which is a scalar quantity) in addition to 
the X-, y- and z-components of the particle motion. 

Many different paths exist by which seismic energy may travel firom the source 1 to a 
receiver 4 in the seismic surveying airangement of Figure 1. A number of paths are 
indicated schematically in Figure 1 . 
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The path 6 shown in Figure 1 is known as the "direct path". Seismic energy that travels 
along the direct path 6 travels fiom the source 1 to a receiver 4 essentially in a straight - 
line without undergoing reflection at any intraf ace. 

Path 7 in Figure 1 is an example of a "water layer multiple path". Seismic energy that 
follows a water layer multiple path propagates wholly within the water column 2, but 
undergoes one or more reflections at the surface of the water column and/or the seafloor 
5 so that the seismic energy passes through the water column more than once. The 
water layer multiple path 7 shown in Figure 1 involves one reflection at the seafloor 5 
and one reflection at the surfece of the water column, but many other water layer 
multiple paths exist 

The path 8 in Figure 1 is an example of a "critical refraction path". Seismic energy that 
foUows the path 8 propagates downwards to tiie seafloor 5, and penetrates into the 
earth's interior 10 (ie into the seabed). The seismic energy continues propagating 
downwardly, untU it reaches a boundary 1 1 between two layers of the earth that have 
different acoustic impedance. The seismic energy undergoes critical refraction, 
propagates along the boundary 11, before eventually being refracted upwards towards . 
the receiver 4. Critical refraction may also occur at the water-seabed interface, and 
downwardly propagating seismic energy that is refracted in this way will propagate 
along the water-seabed interface and wiU then propagate upwardly into the water 
column. 

The path 9 shown in Figure 1 is known as a ''primary reflection path". Seismic energy 
that follows the primary reflection path 9 propagates downwards through the water 
column, is refracted at flie seafloor 5. and propagates downwardly through the earth's 
interior. The seismic energy is reflected at the boundary 11, but is not critically 
refracted and so continues to propagate downwardly into the earth. It eventually 
undergoes reflection at a geological structure 12 that acts as a partial reflector of seismic 
energy, and the reflected seismic energy is, after further refraction as it passes upwardly 
through the boundary 1 1, incident on the receiver 4. The general intent of a seismic • 
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survey is to make use of the seismic energy that foUows fhe primary reflection path in 
order to obtain information about the interior structure of the earth. 

Seismic energy acquired at a receiver may contain iq)wardly and/or do-wnwardly 
propagating seismic energy depending on the location of llie receiver and on the event 
For example seismic energy that travels along the critical refraction path 8 shown in 
Figure 1 will, when it is incident (travelling upwardly) on the water-seabed interface, be 
partly transmitted into the water column and partially reflected back into the seabed 10. 
Thus, a critical refraction event will consist purely of upwardly propagating seismic 
energy above the seafloor 5. but wiU contain botii iq)waidly and downwardly 
propagating seismic energy below Ihe seafloor 5. As another example, seismic energy 
that travels along the direct path 6 shown in Figure 1 will, when incident on the water- 
seab«i interface 5, be partially transmitted into the seabed and partially reflected back 
into the water column. Hence, the direct event will contain both upwardly and 
downwardly propagating seismic energy above the seafloor. but will contain only ' 
downwardly propagating seismic energy below the seafloor. It is therefore often of 
interest to decompose the seismic data acquired at the receiver 4 into an iq)-going 
constitixent and a down-going constituent, above or below the seafloor 5. For example, 
in a 4-C seismic survey it may be of interest to decompose the pressure and tiie vertical 
particle velocity recorded at tiie receiver into their up-going and down-going - 
constitiients above flie seafloor. Various filters that enable decomposition of seismic 
data into i^-going and down-going constituents have been proposed. One example can 
be found in K.M. Schalkwijk et al, "AppUcation of Two-Step Decomposition to Multi- 
Component Ocean-BottomData: Theory and Case Study". J. Seism. Expl. Vol. 8 
pp261-278 (1999), and states that flie down-going and iq>-going constitiients of the 
pressure just above tiie seafloor may be expressed as follows: 
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where P is the pressure acquired at the receiver, is the up-going constitueiit of the 
pressure above the seafloor, is flie down-going constituent of the pressure above the 
seafloor,/is the frequency, itis the horizontal wavenumber, Z is the vertical particle 
velocity component acquired at the receiver.j? is the density of the water, and g is the 
vertical slowness in the water layer. 

As can be seen, the expressions in equation (1) require two of the components of 
seismic data recorded at the receiver to be combined. These filters are an example 
where it is necessary to combine two components of tiie acquired seismic data. It may 
also be necessary to combine two or more components of the acquired seismic data in 
order to decompose the acquired seismic data into p-wave and s-wave (pressure-wave 
and shear-Wave) components, or to remove water level multiple events from the seismic 
data. 

One problem in combining different coinponents of the seismic data acquired at a 
receiver is liiat the different components of the seismic data may not be correctly 
caUbrated against one another. This is particularly the case where the two components 
that are being combined are. as in equation (1), the pressure and the vertical particle 
velocity. There are usuaUy differences in coupling or impulse response between the 
hydrophone used to acquire the pressure and the geophone used to acquire the vertical 
particle velocity. It is necessary to caUbrate the data for these differences before the 
pressure and vertical pdrticle velocity can be combined. This may be done by 
developing a caKbration filter that compensates for the differences in coupling and 
impulse response between the hydrophone and fbe vertical geophone. 

Schalkwijk et al, and others, have suggested that the calibration problem can be 
addressed by assuming that one component of the seismic data has been correctly 
recorded, and caUbrating the other component of the seismic data against the component 
that is assumed to be correctly recorded. In general, it is assumed that the hydrophone 
is weU coupled, so that the pressure recording is taken to be correct The vertical 
component of the particle velocity is then calibrated against the pressure to compensate 
for coupling and impulse response differences between the hydrophone and the vertical 
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geophone. Schalkwijk et al therefore proposed that equation (1) above should be 
modified by applying a caUbration filter to the vertical particle velocity. They proposed 
that the equation given above for the down-going constituent of the pressure above the 
seafloor should be modified to read as follows: 



In equation (2) a(f) represents a frequency-dependent caUbration filter. The remaining 
terms in equation (2) have the same meaning as in equation (1). 

The method proposed by Schalkwijk et al for determining the caUbration filter a(f) is to 
minimise the energy of Ibe down-going pressure constituent above the seafloor for a 
portion of the seismic data that contains only primary reflections. Seismic energy 
travelUng along a primary reflection path is propagating upwardly just above the 
seafloor at tiie receiver position, so that the down-going constituent of the pressure just 
above the seafloor should be zero for data that contains only primary reflections. 
Schalkwijk proposed that the caUbration filter that minimises the energy of flie down- 
going pressure in a window containing only primary reflection events can be found 
using a least squares method. Once the caUbration filter a(f) has been determined in fliis 
way, it is ^pUed to the entire data set for tiie vertical particle velocity. 

The existence of various paths of seismic energy from tiie source to ttie receiver means 
that the data acquired at the receiver in a real seismic survey will contain events 
corresponding to more tiian one possible path. These events will occur at different 
times after tiie actuation of tiie seismic source 1, as different patiis of seismic energy 
have dififerent associated travel times. Figure 2 is a schematic illustiration of seismic 
data tiiat might be acquired at tiie receiver 4. and it shows tiie anq)Utiide of seismic 
energy recorded at tiie receiver 4 as a fimction of tiie time since tiie actuation of flie 
source 1. Figure 2 illustrates a direct event 13, corresponding to the direct patii 6, a 
critical refraction event 14 corresponding to tiie critical reflection patii 8, a primary 
event 15 corresponding to tiie primary reflection patii 9, and a water layer multiple 
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evaQt 16 corresponding to the water layer multiple path 7. (In practice, data acquired at 
a receiver will contain a plurality of primary reflection events from different geological 
structures, a pluraKty of critical refraction events, and a pluraHty of water level multiple 
events arising from different water level multiple paths. Only one event of each type is 
shovm. in Figure 2 for simplicity of explanation.) In order to apply the method of 
Schalkwijk et al to determine the calibration filter, data in a time window that contains 
only the primary event 15, such as the time window A shown in Figure 2, must be 
selected. 

The present invention provides a method of processing multi-component seismic data 
obtained from seismic signals propagating in a medium, the method comprising liie 
steps of: selecting a first portion of the seismic data containing only events arising from 
critical refraction of seismic energy; and determining a first calibration filter from the 
first portion of the seismic data, the first calibration filter being to calibrate a first 
component of the seismic data relative to a second component of the seismic data. 

The method proposed by SchaUcwijk et al has the disadvantage that the time window 
containing only primary reflection events has to be picked manually. The primary 
reflection events are not the first events acquired at the receiver following actuation of 
the source, and so cannot be picked automaticaUy. A further disadvantage is that in 
some cases, for example if the seismic source has a long signature, it may be hard to 
distinguish between the direct arrival and the primary reflection events, so that it may be 
difficult to isolate the correct evcaits. The direct evait contains downwardly 
propagating seismic energy so that use of a time window that inadvertently included the 
direct event would not give correct results for the calibration filter, since the method for 
determining the calibration filter assumes that the selected data contains only iq)-going 
energy. A fiirther problem with the method of Schalkwijk et al is that in shallow water 
the water layer multiple events may arrive at substantially the same time as the primary 
reflection events, and this again makes it difficult to pick a time window that includes 
only "ttie primary reflection events. 
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The present invention makes use of the fact that flie critical refraction events consist 
only of up-going seismic energy just above the seafloor. Thus, selecting a time window 
that contains only one or more critical refraction events makes it possible to detennine 
the caUbratioh filter a(f) by the technique of minimising the energy of the down-going 
pressure just above the seafloor in that time window. 

The method of the invention is particularly advantageous when appUed to long offeet 
data. As is shown in Figure 3, as the offeet (that is/the horizontal distance between the 
source and the receiver) increases, the arrival time of the first critical refraction event 
increases more slowly than does the arrival time of the direct event. For offeets greater 
than Oi the first arrival at the receiver is not the direct event, but is the critical refraction 
event That is. at long offeets the critical refraction event 14 in Figure 2 (two critical 
refraction events are shown in Figure 3) will arrive before the direct event 1 3 and will 
be the first arrival at the receiver. When the invention is appUed to data having a source- 
receiver offeet sufficiently large for the first event acquired at the receiver to be a 
critical refraction event, it is possible to use a time window that covers only the first 
event acquired at the receiver - and tills makes it possible to use an automatic picking 
metiiod to deteraiine tiie time window. If several critical refraction events arrive at the 
receiver before flie direct event, as for far offsets in Figure 3 where two critical 
refraction events arrive before tiie direct event arrives, then all tiiese critical refraction 
events may be included in the time window. 

Figure 4 is a schematic illustiation that coiresponds to Figure 2, but illustrates tiie 
arrival times of tiie events at an offeet tiiat is sufficiently large such tiiat flie first arrival 
is a critical refraction event In tiiis case, flie invention may be appUed by selecting flie 
time vnndow B which includes flie critical refraction event only, and minimising flie * 
energy of flie down-going pressure above flie seafloor in fliis time wmdow. 

A fiirther advantage of tiie invention is tiiat tiie mefliod may be applied to seismic data 
acquired in shallow waters. Altiiough water layer multiple events in seismic data 
acquired in shallow water may coincide witii primary reflection events, tiiey do not 
coincide wifli critical refraction events. Choosing a time window tiiat includes only flie 
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critical refraction event therefore ensures that the time window cannot contain water 
layer multiple evraits. The invention also overcomes the problems that arise when a 
seismic source having a long source signature is used, 

A preferred embodiment of the invention comprises the furftier steps of selecting a 
second portion of tibe seismic data containing only events arising from primary 
reflection of seismic energy and determining a second caUbration filter from the second 
portion of the seismic data, the second calibration filter being to calibrate the first 
component of the seismic data relative to the second component of the seismic data. It 
may comprise the further step of determining a wavenumber-dependent calibration filter 
from the first calibration filter and the second calibration filter. 

A funher problem related to the mefliod proposed by Schalkwijk et al is that the correct 
caUbration filter a(f) may well be dependent on the wavenumber as well as on the 
frequency. The caKbration filter proposed by Schalkwijk, however, is dependent only 
on firequency and, fiirthermore, is derived purely from seismic data at low 
wavenumbers. In an embodiment of the present invention, the filter obtained from the 
critical refraction events is combined with a filter obtained fsova. primary reflection 
events, and a wavenumber-dependent filter is obtained firom the two individual filters. 
The wavenumber-dependent filter may be obtained by, for example, interpolation 
between the filter derived from die critical refiraction eveats and the filter derived fix)m 
the primary reflection events. 

A second aspect of the present invention provides a method of processing multi- 
component seismic data obtained from seismic signals propagating in a medium, the 
metiiod comprising flie step of selecting a first portion of the seismic data corresponding 
to a first wavenumber range; determining a first calibration filter from the first portion 
of the seismic data; selectmg a second portion of the seismic data corresponding to a 
second wavenumber range difEerent from the first wavenumber range; determining a 
second caKbration filter from the second portion of the seismic data; and determining a 
wavenumber-dependent caUbration filter from the first caUbration filter and tiie second 
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calftration filter, the wavanimber-dependent calibration filter being to caHbrate a first 
component of the seismic data relative to a second component of the seismic data. 

A third aspect of the invention provides a method of processing multi-component 
seismic data obtained from seismic signals propagating in a medinm, the method 
comprising the steps of: selecting a first portion of the seismic data in which the first 
arrival contains only upwardly propagating seismic energy above the seafloor; and 
detennining a first caUbration filter fi-om the first portion of the seismic data, the first 
caUbration filter being to calibrate a first component of the seismic data relative to a 
second component of the seismic data. 

The invention may be applied to any event that is the first arrival and that contains only 
iipgoing energy above the seafloor. For example, at fer offsets tiie first arrival may be 
an event which is not a critical refraction event but v^Mch nevertheless contains only 
upgoing energy above the seafloor - such as, for exanq)le, a wave that was trapped in a 
tiiin subsurfece layer of the seabed - and tiie inventionmay be appUed to such events. 

The invention may fiorther comprise the step of caHbrating the first component of the 
seismic data using the first caUbration filter or using the wavenumber-depettdent 
caUbration filter. 

A fourth aspect of the present invention provides a method of seismic surveying 
comprising the steps of: actuating a source of seismic energy; acquiring seismic data at 
a receiver spatiaUy separated from the source; and processing the seismic data by a 
method as defined above. 

A fifth aspect of tiie present invention provides an apparatus for processing multi- 
component seismic data to determine a caUbration filter to caUbrate a first component of 
tiie seismic data relative to a second component of tiie seismic data, the apparatiis 
comprising: means for selecting a first portion of tiie seismic data containing only 
events arising from critical refraction of seismic energy; and means for detennining a 
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first calibration filter from the first portion of the seismic data. The apparatus may 
comprise a programmable data processor. 

A sixfti aspect of tiie invention provides an apparatus for processing multi-component 
seismic data to detemiine a caHbration filter to caUbrate a first component of tiie seismic 
data relative to a second conq)onent of tiie seismic data, the ^paratus comprising: 
means for selecting a first portion of the seismic data in which the first arrival contains 
only upwardly propagating seismic energy above the seafloor; and means for 
deteraiining a first caUbration filter from the first portion of tiie seismic data. 

A seventh aspect of the invention provides an apparatus for processing multi-component 
seismic data to detemiine a caHbration filter to caHbrate a first component of the seismic 
data relative to a second conq>onent of the seismic data, the apparatus comprising: 
means for selecting a first portion of the seismic data corresponding to a first 
wavenumber range; means for deteraiining a first caUbration filter fix>m the first portion 
of the seismic data; means for selecting a second portion of the seismic data 
corresponding to a second wavenumber range different fix>m tiie first wavenumber 
range; means for deteraiining a second calibration filter &om the second portion of the 
seismic data; and means for determining a wavenumber-dependent calibration filter 
from the first calibration filter and tiie second calibration filter. 

The sqpparatus may conq>rise a programmable data processor. 

An eighth aspect of the present invraition provides a storage medium containing a 
program for an ^^aratus as defined above. 

The invention also provides a method of deteraiining a first calibration filter for 
calibrating a first component of multi-component seismic data relative to a second 
component of tiie seismic data, flie metiiod comprising tiie steps of: selecting a first 
portion of the seismic data containing only events arising from critical refiaction of 
seismic energy; and determining the first calibration filter froin the first portion of the 
seismic data. 
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The mvention also provides a method of detennining a wavenmnber-dependent 
calibration filter for calibrating a first conq)onent of multi-component seismic data 
relative to a second component of the seismic data, the method comprising the steps of: 
selecting a first portion of the seismic data corresponding to a fiaSt waveamnber range; 
detennining a first caUbration filter &om the first portion of the seismic data; selecting a 
second portion of the seismic data coiresponding to a second wavenumber range 
different firom the first wavenumber range; detennining a second cahbration filter &om 
the second portion of the seismic data; and determining a wavenumber-dependent 
calibration filter firom the first calibration filter and the second caUbration filter. 

The invention also provides amethod of deteraiining a first calibration filta: for 
calibrating a first component of multi-conq)onent seismic data relative to a second 
component of the seismic data, flie method comprising the steps of: selecting a first 
portion of the sasmic data in which the first arrival contains only Tq>wardly propagating 
seismic energy above the seafloor; and deteraiining a first calibration filter &om the first 
portion of the seismic data. 

Prefenred embodiments of the present invention will now be described by way of 
illustrative example with reference to the accompanying figures in which: 

Figure 1 is a schematic illustration of a seismic survey; 

Figure 2 is a schematic illustration of tiie seismic energy acquired at a receiver in the 
sasmic survey of Figure 1; 

Figure 3 is a schematic illustration of the variation of arrival time of seismic energy as a 
fimction of offset between flie source and the receiver; 



Figure 4 is a schematic illustration of the seismic energy acquired in tiie seismic survey 
of Figure 1 at long offsets, illustrating a method of the present invention; 
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Figure 5 is a schematic illustration of pressure recorded at a receiver in the seismic 
surveying arrangement of Figure 1; 



Figures 6 and 7 illustrate the up-going and down-going constituents of pressure above 
the seafloor obtained from the pressure data shown in Figure 5 according to a prior art 
approach; 

Figures 8 and 9 illustrate the up-going and down-going pressure constituents above the 
seafloor obtained from the pressure data of Figure 5 according to a method of the 
present invention; 

Figure 10 is a schematic block flow diagram of a method of the present invention; and 

Figure 1 1 is a block schematic diagram of an apparatus according to the present 
invention. 

Figure 5 illustrates typical pressure data recorded at a 4-C receiver in a seismic survey 
such as the survey shown in Figure 1. The x-axis in Figure 5 indicates the offeet 
between the source and the receiver, and the y-axis indicates the time after actuation of 
the seismic source. The data are common receiver data and were acqufred using a 
single receiver and a linear array of sources deployed with a spacing of 50m between 
each pair of adjacent sources. Each frace represents the pressure acquired at the 
receiver when one source is actuated, with Ihe an^htude of the acquired pressure being 
in the x-direction. 

It should be noted that different receivers in an anay may well have different coupling, 
different instruments responses etc, even if all the receivers are nominally identical to 
one another. The calibration filter required for data acquired at one receiver m a 
receiver array is therefore likely to be different from tiie caHbration filter required for 
data acquired at another receiver m the array. The invention is therefore preferably 
applied to common receiver gathers and a separate calibration filter is determined for 
each common receiver gather. 
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The pressure data shown in Figure 5 contains a large nmnber of seismic events. The 
event labelled 13 is the direct wave, and it wiU be seen that this is the first arrival for 
o^ets having a magnitude of up to proximately 1,000m. The event labelled 14 is a 
critical refraction event, and it will be seen that this is the first arrival for offsets having 
a magnitude significantly greater than 1,000m. 

Figures 6 and 7 illustrate the up-going constituent above the seafloor (Figure 6) and the 
down-going constituent above the seafloor (Figure 7) of the pressure shown in Figure 5 
obtained using the filters given in equation (1) above. That is, the up-going and down- 
going constituent shown in Figures 6 and 7 were obtained on the assumption that the 
pressure data and the vertical particle velocity data (not illustrated) were correctly 
caKbrated to one another. Inspection of Figures 6 and 7 shows that this assumption is 
incorrect In particular, the critical refi:action event 14 contains only up-going energy 
above the seafloor and so should appear only in the up-going pressure constituent and 
should not pear in the down-gomg pressure constituent It will, however, be seen lhat 
the up-going critical refraction event has leaked through into the down-gomg pressure 
constituent shown in Figure 7, and this indicates that the caHbration is unsatisfectory. 

According to the present invention, a calibration filter is determined from the critical 
refraction event 14. As noted above, for traces acquired at a source-receiver offset 
having a magnitude well above lOOOm, the critical reflation event is the first event 
acquired at the receiver, and is well-separated from the subsequent event. It is therefore 
possible for such traces to define a time-ofiFset window that includes only the first 
critical refraction event, and so includes only up-going energy. 

One suitable time-offset window of data is illustrated in Figure 5, as the region C. It 
will be seen that this region includes traces acquired at an oflfeet of between -3000m to 
approximately -2100m. For each trace in this offset range the region C defines a time 
window that includes only the first refraction event (which is the first arrival in each of 
the selected traces). It will be noted that tbe centre point of the time window for a 
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paiticular trace is not constant betweesn traces but increases with increasing magnitude 
of ofEset. 

The caHbration filter for the vertical velocity component is then calculated on the 
assumption that the energy in the selected portion C of the data should contain only up- 
going energy. The calibration filter may be determined in any suitable way. In 
particular, the cahTjration filter a(f) may be determined by finding the caHbration filter 
that minimises the energy of the down-going pressure constituent using a least squares 
process, as in the method of Schalkwijk et al. Once the appropriate caHbration filter a(f) 
has been determined, revised filters for determining the up-going and down-going 
constituents of the pressure above the seafloor can be determined using equation (2), or 
in general the filter a(f) may be i^pHed to the entire gather of the vertical conq)onent 
data, and the caHbrated vertical component data can then be used as an input to any 
process requiring a combination of the vertical component with any other seismic 
corD|)onents. 

Another suitable portion of data exists in the corresponding region for positive offeets in 
the range 2100 to 3000m. One possible implementation of flie method would be to use 
both these regions, by defining a second region, analogous to the region C in Figure 5, 
for offsets in tiie range +2100m to +3000m and determining a second caHbration filter. 
The two filters determined firom flie two windows may flien be averaged. This will 
however not be possible for all data sets since a receiver gather does not necessarily 
have the same amount of positive and negative oflfeets, and hence a region with a clearly 
separated critical refraction event may be present only for eitiier positive or negative 
of&ets. 

Figures 8 and 9 iUustiate the results of decomposmg the pressure data of Figure 5 into 
its up-going and down-going constituents above the seafloor using filters of the type 
given in equation (2) above, and with the calibration filtOT a(f) determined fix>m the 
seismic data in the region C of Figure 5. It wiU be noted that the critical refiiaction 
event 14 appears predominantly in the up-going pressure constituent of Figure 8, and is 
ahnost completely absent from the down-going pressure constituent of Figure 9. This 
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iUustrates that flie decomposition of Figures 8 and 9 is significantly more accurate than 
the decomposition of Figures 6 and 7, since Ihe critical refraction event is expected to 
occur only in the up-going pressure constituent 

It will also be noted that the primary reflection event is stronger in the up-going 
pressure constituent of Figure 8 than in the up-going pressure constituent of Figure 6. 
This suggests that the cahbration filter found from the critical refraction event at long 
ofEsets is also applicable at low offsets. 

Figure 10 is a block flow diagram illustrating one embodiment of the method of the 
present inv^tion. 

MtiaUy, at step 17, seismic data is acquired. This may be, for example, acquired in a 
survey of the type shown in Figure 1 . 

The invention may alternatively be appUed to.pre-existmg seismic data. Step 17 may 
therefore be replaced by the alternative step 18 of retrieving pre-existing seismic data 
from storage. 

At step 19 a suitable offset range is selected. In the example described above with 
reference to Figure 5, step 19 consists of selecting the offeet range from -3,000m to - 
2100m. 

At step 20, the first arrival of seismic energy for each trace in the selected offset range is 
determined (this may be thought of a selecting a time window for each trace, and so 
defining an of&et-time window). Assummg that flie'offeet range has been selected 
correctly in step 19, the first arrival in each trace in the selected offeet range will be a 
critical refraction event such as the event 14. Since the event is the first event in each 
trace, step 20 may be carried out using an automatic picking method, although it may 
alternatively be perfonned manually. 
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At Step 21, a calibration filter is detennined that is the best fit to the data in the selected 
offset range and time window. This is done by calculating title down-going pressure 
constituent above the seafloor from the pressure and vertical particle velocity recorded 
at the receiver using equation (2), and finding the cahbration filter that minimises the 
energy in the down-going constituent of the pressure. 

At step 22 the filter a0 is appKed to all the desired traces of the vertical component of 
the seismic data acquired at step 17 or retrieved fix>m storage at step 18. 

At step 23 the caUbrated vertical component data is then used as input into any process 
requiring a combination of several seismic components. For example, filters for 
detennining the up-going and down-going constituents of tihe pressure above the 
seafloor may be determined, usmg equation (2) and flie calibration filter determined at 
ste^21. 

If desired, steps 22 and 23 may be omitted. In tiiis case the cahbration filter determined 
at step 21 may be output for display or stored for subsequrat use. 

hi an altemative embodiment of the invention, a wavenumber-dependent filter is 
determined by combining tiie approach of Schalkwijk et al with the present invention, 
hi tiiis embodiment, a cahbration filter is determined firom critical refiraction events 
occurring at long ofifsets, as described above with reference to steps 17 to 21 of Figure 
10. A second cahbration filter is tiien determmed firom traces m which it is possible to 
define a time-offeet wmdow that contams only primary reflection events. A cahbration 
filter is determined for these traces m the manner described by Schalkwijk et al. A 
suitable region of data for obtaming this filt«- is indicated on Figure 5 as region D. 

The cahbration filter determmed firom critical refiaction events at long offeet and the 
cahbration filter determmed firom primary reflection events at low offeet are then 
combmed to produce a wavenumber-dependent cahbration filter. The filters may be 
combined using an mterpolation technique to determme the filter to be apphed at a 
given offeet. 
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In this embodiment, step 22 of Figure 10 is replaced by the step of calibrating tiie 
vertical component using the wavenumber-dependent calibration filter. Alternatively 
steps 22 and 23 may be omitted, and the wavenumber-dependent caKbration filter can 
be oulput or stored for fixture use. 

An alternative way to obtain a wavenumber-dependent caHbration filter is to conq)ute a 
caUbration filter for each separate trace in the offset range selected at step 19. In this 
■ alternative embodiment steps 20 and 21 are performed on each trace (or on a pluraHty of 
selected traces) in the offset range selected at step 19 so that calibration filters are 
determined for several different wavenumbers. Alternatively, the traces in the offset 
range selected at step 19 can be groined, and a calibration filter can be determined for 
each groiq) of traces, for example usmg a least squares method. Again, this results in 
calibration filters for several different wavenumbers. 

Once calibration filters have been obtained for several different wavenumbers, it is 
possible to interpolate between and/or extrapolate from these caKbration filters to obtain 
a wavenumber-dependent caHbration filter. This method would, however, only work 
well for a time-offset wmdow containing only primary reflections (i.e., the window D in 
Figure 5), since the vertical slowness is constant for the refiracted event. The 
wavenumber-dependent caUbration filter may be used immediately, or may be output or 
stored for future use. 

A fiirflier alternative method is to define time-of&et windows around several refiraction 
events of different vertical slownesses, and determine aplurality of caUbration filters 
(one caUbration filter can be obtained &om data in each window). A wavenumber- 
dependent caUbration filter can be obtained by interpolation between and/or 
extrapolation from these caUbration filters. If desired, one or more caUbration filters 
determined from a time-of£set window containing only primary reflections can also be 
used in the interpolation and/or extrapolation. The wavenumber-dependent caUbration 
filter again may be used immediately, or may be output or stored for fixture use. 
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It wiU be noted in Figure 5 that the primary reflection event is obscured by other events 
at long offeets. It will therefore be extremely difficult to compute a reUable calibration 
filter at long offeets using the method of Schalkwijk et al, owing to the difficulty of 
detennining a time window that contains only primary reflection events. Furthennore, 
even if a time window that contained only primary reflection events could be 
detennined for the long offset traces in Figure 5, fliis could only be done by a manual 
picking mediod and could not be automated. 

The invention has been described above with reference to a calibration filter that 
calibrates the vertical particle motion with regard to the pressure, on the assumption that 
the pressure has been accurately recorded. The invention is not limited to this, however, 
and in principle could be used to determine a calibration filter that calibrates the 
pressure with regard to the vertical particle motion, on the assumption that the vertical 
particle motion has been accurately recorded. 

Figure 1 1 is a schematic block diagram of an ^paratos 34 according to the present 
invention. The apparatus is able to cany out a method according to the present 
invention. 

The apparatus 34 comprises a programmable data processor 27 with a program memory 
28, for instance in the form of a read only memory ROM, storing a program for 
controlling the data processor 27 to process seismic data by a method of the invention. 
The ^paratus fiirther comprises non-volatile read/write memory 29 for storing, for 
example, any data which must be retained in the absence of power supply. A ' Vorking" 
or "scratclq>ad" memory for the data processor is provided by a random access memory 
(RAM) 30. An input device 31 is provided, for instance for receiving user commands 
and data. An output device 32 is provided, for instance for displaying information 
relating to the progress and result of tiie method. The ou^ut device may be, for 
example, a printrar, a visual display unit or an ou^ut memory. 

Seismic data for processing may be supplied via the input device 31 or may optionally 
be provided by a machine-readable store 33. 
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The program for operating the ^axatus and for performing a method as described 
hereinbefore is stored in the program memory 28, which may be embodied as a semi- 
conductor memory, for instance of the well-known ROM type. However, the program 
may be stored in any other suitable storage medium, such as magnetic data carrier 28a 
(such as a "floppy disc") or CD-ROM 28b. 
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CLAIMS: 

1 . A method of processing multi-coniponeat seismic data obtained from seismic 
signals propagating in a medium, the method comprising the steps of: selecting a first 
portion of the seismic data containing only events arising from critical refraction of 
seismic energy; and determining a first calibration filter from the first portion of the 
seismic data, the first calibration filter being to calibrate a first component of the 
seismic data relative to a second component of the seismic data. 

2. A method as claimed in claim 1 wherem the first portion of the seismic data is 
data acquired with a long source-receiver offset. 

3. A method as claimed in claim 1 or 2 wherein the fijrst component is the vertical 
component of particle motion and the second component is pressure, 

4. A method as claimed in claun 1 or 2 wherem the first component is pressure and 
the second component is the vertical component of particle motion. 

5. A method as claimed in claim 1, 2, 3 or 4 wherein the step of detemaining &e 
first calibration filter comprises minimising the energy immediately above the seafloor 
of the downgoing constituent of the second component for the selected portion of the 
seismic data. 

6. A method as claimed in any preceding claim and comprising the finiher steps of 
selectmg a second portion of the seismic data containing only events arising from 
primary reflection of seismic energy and determining a second calibration filter from the 
second portion of the seismic data, the second calibration filter being to calibrate the 
first component of the seismic data relative to the second component of the seismic 
data. 
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7. A method as claimed in claim 6 and con^rising the further step of determining a 
wavenumber-dependent calibration filter firom the first caHbration filter and the second 
calibration filter. 

8. A method of processing multi-component seismic data obtained &om seismic 
signals propagating in a medinm, the method comprising Hie steps of: selecting a first 
portion of liie seismic data corresponding to a first wavenumber range; detennining a 
first caHbration filter firom the first portion of the seismic data; selecting a second 
portion of the seismic data corresponding to a second wavenmnber range different firom 
the first wavenumber range; determining a second caUbration filter from the second 
portion of the seismic data; and detennining a wavenumber-dependent caHbration filter 
from the first caHbration filter and the second caHbration filter, the wavenumber- 
dependent caHbration filter being to caHbrate a first component of the seismic data 
relative to a second component of the seismic data. 

9. A method as claimed in claim 8 wherein the first wavenumber range 
corresponds to seismic data containing substantiaUy only critical refiraction events and 
the second wavenumber range corresponds to seismic data containing substantiaHy only 
primary reflection events. 

10. A method of processing multi-component seismic data obtained from seismic 
signals propagating in a medium, the method comprising the steps of: selecting a first 
portion of the seismic data in which the first arrival contains only upwardly propagating 
seismic energy above the seafloor; and deteraiimng a first caHbration filter from the first 
portion of the seismic data, the first caHbration filter being to caHbrate a first component 
of the seismic data relative to a second component of the seismic data. 



11. A method as claimed in any preceding claim and comprising the further step 
caHbrating the first component of the seismic data using the first caHbration filter. 
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12. A method as claimed in claim 7, 8 or 9 and comprising the further step of 
calibrating the first component of the seismic data using the wavenumber-dependeut 
calibration jQlter. 

13. A method of seismic surveying comprising the steps of: actuating a source of 
seismic energy; acquiring seismic data at a receiver spatially separated firom the source; 
and processing the seismic data by a method as defined in any of claims 1 to 12. 

14. An apparatus for processing multi-component seismic data to determine a 
calibration filter for calibrating a first component of the seismic data relative to a second 
component of the seismic data, the apparatus comprising: means for selecting a first 
portion of the seismic data containing only events arising &om critical refi-action of 
seismic energy; and means for determining a first calibration filter firom the first portion 
of the seismic data. 

15. An jQ)paratus for processing multi-component seismic data to determine a 
calibration filter for caHbrating a first component of the seismic data relative to a second 
component of the seismic data, flie apparatus comprising: means for selecting a first 
portion of the seismic data in which the first arrival contains only inwardly propagating 
seismic energy above the seafloor; and means for determining a first calibration filter 
firom the first portion of the seismic data. 

16. An apparatus for processing multi-component seismic data to determine a 
calibration filter for calibrating a first component of the seismic data relative to a second 
component of the seismic data, flie ^aratus comprising: means for selecting a first 
portion of the seismic data corresponding to a first wavenumber range; means for 
detemiining a first calibration filter fimm the first portion of flie seismic data; means for 
selecting a second portion of the seismic data corresponding to a second wavenumber 
range different fix>m liie first wavenumber range; means for dete rmining a second 
caUbration filter firam the second portion of the seismic data; and means for determining 
a wavenumber-dependent calibration filter firom the first caUbration filter and the 
second caUbration filter. 
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17. An ^aratus as claimed in claim 14 or 15 and fiirfher comprising means for 
calibrating the first component of the seismic data using the first calibration filter. 

18. An ^paratus as claimed in claim 16 and fbrfher conq>rising means for 
cahT)rating the first component of the seismic data using the wavenumber-dependent 
calibration filter. 

19. An apparatus as claimed in any of claims 14 to 18 and comprising a 
programmable data processor. 

20. A storage medium containing a program for an apparatus as defined in claim 
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ABSTRACT 

A method of and apparatus for processing seismic data 



A method detemuning a calibration filter to calibrate a first component of multi- 
component seismic data relative to a second component of the seismic data comprises 
detemiining the calibration filter from a portion of the seismic data that contains only 
events arising from critical refraction of seismic energy. The method is particularly 
suitable for long-offset data, since the first arrival will be a critical reflation event and 
an automatic picking method may be used. 

The present invention also provides a wavenumber-dependent calibration filter that is 
obtained from a calibration filter obtained from data in one offset range and another 
calibration filter obtained from data in another offset range. 
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Figure 5: Raw pressure data with two possible minimlzatioii windows; one containing the critical 
refraction event (C), and one containing the primary reflections (D). 
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Figure 4: Upgoing pressure when no calibration filter is applied. 
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Figured: Downgoing pressure when no calibration Olter is applied* 
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Figured Upgoing pressure when calibration Alter is applied. 
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Figured: Downgoii^ pressure when calibration filter is applied. 
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